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^ ' Abstract 

o ; 

■ We analyse subleading electroweak effects in the top anti-top resonance produc- 

tion region in e"'"e~ collisions which arise due to the decay of the top and anti-top 
quarks into the W~^W~bb final state. These are NLO corrections adopting the non- 
^ ' relativistic power counting ~ ~ \/oew- In contrast to the QCD corrections 

which have been calculated (almost) up to NNNLO, the parametrically larger NLO 
electroweak contributions have not been completely known so far, but are manda- 
tory for the required accuracy at a future linear collider. The missing parts of these 
NLO contributions arise from matching coefficients of non-resonant production- 
decay operators in unstable-particle effective theory which correspond to off-shell 
top production and decay and other non-resonant irreducible background processes 
to production. We consider the total cross section of the e^e~ — >• W~^W~bb pro- 
cess and additionally implement cuts on the invariant masses of the W~^b and W~b 
pairs. 



1 Introduction 



The top quark is currently known from direct production at the Fermilab Tevatron 
to weigh nit = 173.1 ± 0.6 (stat.) ± 1.1 (syst.) GeV [T], and an increase in precision to 
about 1 GeV is expected soon from the Large Hadron Colhder (LHC). From a thresh- 
old scan of the e'^e~ — )■ tt cross section at the planned International Linear Collider 
(ILC), however, a precision of 30MeV can be achieved experimentally |2]. Aside from 
determining a fundamental parameter of the Standard Model, accurate top-mass mea- 
surements constrain the quantum fluctuations from non-standard interactions in elec- 
troweak precision measurements. Other characteristics of the top quark such as its 
width and Yukawa coupling provide information about its coupling to other particles 
and the mechanism of electroweak symmetry breaking. For these reasons (and perhaps 
also because it represents the cleanest non-relativistic system bound by the colour force) 
top-quark pair production near threshold in e'^e~ annihilation has been thoroughly in- 
vestigated following the non-relativistic approach of [3K1I5] . which treats the leading 
colour- Coulomb force exactly to all orders in perturbation theory. In this framework, 
where the strong coupling as is of the same order as v, the small relative velocity of the 
top and anti-top, next-to-next-to-leading order (NNLO) corrections have been available 
for some time PI7|8|9|T0|11|12] . next-to-leading and some higher-order logarithms of v 
have been summed to all orders P^14|15|[TB] . and the third-order (NNNLO) cross sec- 



tion is now known almost completely [T7|18|I19] . which requires input from three-loop 
matching coefficients [20], potentials [2T|22f23f24f25] . and third-order S-wave energy lev- 
els and residues [T7|26f27|l28] . The full NNNLO result should finally clarify the question 
whether the QCD corrections can be calculated with the required precision. 

Here we focus on a different issue, non-resonant production of the physical final state 
after top decay, which though known for some time (see, e.g., the discussion in [29] ) 
has been left aside up to now. The top quark is unstable with a significant width Ft of 
about 1.5 GeV due to the electroweak interaction. The width is essential in threshold 
production, since it prevents the top and anti-top from forming a bound state [30] and 
causes a broad resonance structure in the energy dependence of the cross section on top 
of the increase due to the opening-up of the two-particle phase space. However, once 
the top width is included, due to top decay, the physical final state is W~^W~bb - at 
least if we neglect the decay of top into strange and down quarks, as justified by Vib ~ 1, 
and consider W bosons as stable, which we may assume for the purpose of discussion 
in the present paper. The W^W~bb final state can be reached through non-resonant 
top production or background processes containing no or only single top quarks. Both 
effects are not included in the non-relativistic treatment. Adopting a standard counting 
scheme where oew ~ af, we find that the leading non- resonant and off-shell effects are 
next-to-leading order (NLO) for the total cross section, since there is an additional power 
of oew but no phase-space suppression, hence the relative correction is oew/^ ~ C(s- We 
note that some of the higher-order NNLO corrections due to the finite top width have 
already been calculated in [21]. 

In this paper we calculate the complete NLO "electroweak" contributions to the 
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e^e~ — > W^W~hh process in the ti resonance region, for the total cross section as well 
as including invariant-mass cuts on the Wh pairs. The calculation is performed with 
unstable-particle effective field theory [32|33] . which provides the framework for consis- 
tently including resonant and non-resonant effects while maintaining an expansion in the 
small parameters of the problem. There are many similarities with IV-pair production 
considered in [M1I35] . though the top-quark case is technically more complicated due 
to the presence of a massive particle (the W boson) in top decay. We note that an 
alternative approach has been developed in parallel [36f37] that includes the effects of 
invariant-mass cuts on the Wh pairs entirely through calculations in non-relativistic ef- 
fective theory. This works if the invariant-mass cuts around rrit are neither very loose nor 
very tight, and provided that the non- resonant background processes are small, which 
must be checked by other means. The present approach removes both these restrictions 
and can be used in particular for the total cross section. We shall see that when both 
approaches are applicable, they agree well for top-quark production. 

There is an interesting conceptual issue concerning the pure QCD calculation of 
the tt cross section. It exhibits an uncanceled ultraviolet divergence (here regularized 
dimensionally) 

att oc oc — ^ — (1) 

at NNLO, which arises from the logarithmic overall divergence in the two-loop non- 
relativistic correlation function, whose imaginary part gives the cross section. The over- 
all divergence is polynomial in the non-relativistic energy E of the top quarks, but 
contributes to the cross section, since the correlation function is evaluated at complex 
values E ^ E + iVt- As discussed in [19] this divergence cancels with an infrared diver- 
gence that appears in the non-resonant term in unstable-particle effective theory from 
the diagram corresponding to off-shell top-quark decay. The lesson from this is that the 
pure QCD result alone that is usually shown in the literature is inconsistent theoretically 
and must be embedded in the systematic calculation of the e+e" — )■ W~^W~bb process. 
In the present paper, however, we are concerned with NLO accuracy, where no explicit 
divergence arises in dimensional regularization. But the problem at NLO is in fact worse. 
Inspection of the one-loop non-relativistic correlation function shows that it is linearly 
ultraviolet divergent, or more generally exhibits a cut-off dependence of the NLO or- 
der, just as the non-resonant contribution that we calculate below has a linear infrared 
divergence. Both require regularization and as usual dimensional regularization yields 
finite (but regularization-dependent) values in four dimensions for linearly divergent in- 
tegrals. The conclusion is that the non-resonant NLO contribution that we consider here 
is mandatory to obtain a regularization-independent result for top anti-top production, 
more precisely W~^W~bb production in the tt threshold region, at NLO. 

The outline of the paper is as follows. In Section [2] we briefly review unstable- 
particle effective theory. We list the terms that contribute at NLO to the e~^e~ — )■ 
W~^W~bb process, when counting as ~ a/oew ~ 'V, and explain that the non-resonant 
electroweak contribution is obtained from the e'^e~ — )■ t W~b amplitude and its charge- 
conjugate with the top-quark width set to zero. The computation is described and 
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the results are summarized in Section [3] and Appendix |Al In Section H] we discuss the 
implementation of invariant-mass cuts in the effective-theory calculation and validate 
our result by comparing it to the Born cross section for the e+e~ — )■ W~^W^bb process 
computed numerically with MadGraph/MadEvent. Finally, in Section [5] we combine 
the new non-resonant contributions with the resummed leading-order QCD result and 
the electromagnetic correction due to the electromagnetic Coulomb potential, which 
constitutes another electroweak NLO effect, and discuss the top anti-top line-shape as 
seen in the energy-dependence of the physical W~^W~bb final state compared to the pure 
QCD calculation of ti production. We conclude in Section [61 



2 Method of calculation 

2.1 Unstable-particle effective theory for pair production near 
threshold 

The cross section for the e~^e~ — ?■ W~^W~bb process (inclusive but also in the case where 
we allow for cuts in the invariant masses of the W~^b and W~b subsystems) is obtained 
from the W~^bW~b cuts of the e'^e~ forward-scattering amplitude. In the energy region 
a/s !^ 2mt close to the top anti-top production threshold, the amplitude is dominated by 
the production of resonant top quarks with small virtuality. This allows us to integrate 
out hard modes with scale rrit and represent the forward-scattering amplitude as the sum 
of two terms [32][33l. 




+ J2d::\e-e^\^Ot:^me-e^). (2) 

k 

In the subsequent discussion of top-quark pair production near threshold we follow closely 
the formalism described in [31] for four-fermion production in e~^e~ collisions in the 
energy region of the H^-pair production threshold, and refer to that paper for further 
details on the method. 

The matrix elements in ([2]) are evaluated in the "low-energy" effective theory, which 
includes elements of soft-coUinear and non-relativistic effective theory. The first term 
on the right-hand side of ([2]) describes the production of a resonant tt pair in terms 
of production (decay) operators Op\x) {Op''^\x)) with short- distance coefficients C^'^\ 
The second term accounts for the remaining non-resonant contributions. Note that there 
is no separate term for the production of one resonant and one off-shell top [MEZ], since 
such configurations are not sensitive to the particular low-energy dynamics that develops 
at the pair-production threshold. They are effectively short-distance and included in the 
coefficient functions C^^^ of non-resonant production-decay operators 0^^^{'S). 



3 



The standard pure QCD calculation of top-pair production near threshold is entirely 
contained in the resonant term. In this paper we are mainly concerned with the calcu- 

(k) 

lation of the leading contribution to Cl^\ which represents an NLO correction to the 
forward-scattering amplitude A, since the first term in ([2]) is of order a'^yijV, while the 
second is 0{a^-s^). Nevertheless, for comparison and to discuss other NLO non-QCD 
effects, we briefly summarize the expressions relevant to the leading-order resonant term. 
There are only two production operators at this order. 



p 



^C2 



p 



(3) 



where ipt {xt) denotes the non-relativistic top (anti-top) field, a collinear electron field 
with large, light-like momentum in a direction labelled by c, and 7^ = —7*. We drop the 
Wilson lines required to make the operators gauge- invariant, since they do not contribute 
to our calculation. The tree-level coefficient functions of these operators are computed 
from the tree-level amplitude of the process e^e~ — )■ ti at production threshold i/i = 2mt 
and are given by (keeping the full s-dependence in the photon and Z propagators) 



Op 



Ana 
—Ana 



QtQe 

S S 

s-Mr 



VeVt 



(4) 



where a denotes the electromagnetic coupling, Q j the electric charge of fermion species 
/ in units of the positron charge (Qt = 2/3, Qe = —1), and Mz the Z-boson mass. The 
vector and axial-vector coupling of the fermion / to the Z-boson are given by 



n - 2Qfsl 



T 



^23 11)0^1 



(5) 



where Sy^ (c^) is the sine (cosine) of the weak mixing angle and T/ the third component 
of the weak isospin of the fermion. In terms of these quantities the leading-order cross 
section is given by 



a 



(0) 



-Im^(°) 

s 



where E = ^/s — 2mj and S 



L/p -I- Up 



E + iTu N, 



2iVJmG'[?^(0,0;^), (6) 
3 equals the number of colours, and 



G[?'*(0, 0; £^) denotes the MS-renormalized zero-distance Coulomb Green function [38] 



4°^ (0,0;^) 




(7) 
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Here 7^; is the Euler-Mascheroni constant, /i the scale introduced in dimensional regu- 
larization, and '0 the Euler psi-function. (The dependence on n and the regularization 
drops out in ImG[?^(0, 0; 

The leading contribution from non-resonant production-decay operators OIJ to ([2j) 
arises from four-electron operators of the form 

O'^J = eciTiCca ecaTaCci, (8) 

where Fi, r2 are Dirac matrices. The calculation of the short- distance coefficients C^^^^ 
is performed in standard fixed-order perturbation theory in the full electroweak theory. 
In particular, the top propagator is the free propagator not including the top width, 
which ensures that the amplitude depends only on the short- distance scales. Self-energy 
insertions are treated perturbatively and are not resummed into the propagators. The 
leading contribution to the forward-scattering amplitude arises from the one-loop dia- 
grams with two top propagators. Just as in the case of VT-boson pair production [M] 
the imaginary part of these diagrams vanishes in dimensional regularization to all orders 
in the expansion in E = y/s — 2mt in the hard region. Thus the leading imaginary 
parts of C4g^ arise from two-loop diagrams of order shown in Figure [1] below. The 
corresponding contribution to the cross section is 

(Tnon-res = ^ 5Z [^4e^] {e' e+ \^0^£ . (9) 

k 

It is convenient not to calculate the separately, but directly the sum ([9]) when one 
considers the unpolarized cross section. Technically, this simply amounts to the calcula- 
tion of the spin-averaged tree-level processes e~^e~ — )■ tW~b and e+e~ — )■ tW~^b with no 
width supplied to the intermediate top-quark propagators. Instead, the divergence from 
the top-quark propagators going on-shell is regularized dimensionally. 



2.2 NLO corrections related to the top-quark instability and 
electroweak interactions 

In the following we list the NLO corrections to the e~^e~ — )■ W^W~bb process in the 
top anti-top resonance region ^/s ~ 2mt. We do not discuss in this paper the pure 
QCD NLO corrections consisting of the one-loop QCD correction to the 7*(Z)tt vertex 
and the one-loop correction to the QCD Coulomb potential, since they are standard 
and included in the pure QCD calculations in the literature. Instead we focus on NLO 
corrections which arise once we consider the W~^W~bb final state (as required by the top 
instability), and other electroweak effects: 

• The single insertion of the (leading-order) electromagnetic Coulomb interaction 
instead of the colour Coulomb interaction is an NLO effect, since a/ag ~ a^. 
This is trivially accounted for, including some higher-order terms, by replacing 
asCp — asCp + (yQt in ©• This effect has already been studied in [T6|l37] . 
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• The one- loop QCD correction to the on-shell top width [39] is also trivially included 
by using the one- loop corrected expression for Ft in the variable S = E + iTt in ([7]). 
We do not discuss this effect further in this paper, however, since in the pure QCD 
calculation the top width is considered as an adjustable input parameter, and thus 
may be thought of as representing the one-loop corrected expression. 

• Corrections coming from gluon exchange involving the bottom quarks in the final 
state (QCD interference effects) appear at NLO but vanish at NLO for the total 
cross section [IQ] . We show in Section |4] that they are parametrically suppressed 
in the presence of invariant-mass cuts, provided the cuts are loose in the sense 
defined below (see also [37]). Since we focus on loose cuts in the final results, no 
calculation of interference effects is required. 

• Non-resonant contributions start at NLO, as discussed above. They arise from the 
hard contributions to the Born cross section for e~^e~ — )■ W~^W~bb, which deter- 
mine the matching coefficients of the four-electron pro duct ion- decay operators. We 
calculate this contribution in Section [21 

We should note that any three-particle final state produced through the electroweak 
interaction in e~^e~ annihilation contributes at NLO (namely 0{a^-y^)) provided the 
W^W~bb final state can be reached through the resonant decay of one of the three 
particles with a significant branching fraction of the decaying resonance. The obvious 
relevant cases are t W~b and t W~^b considered in Section[3l others are e~^e~ — )■ W~^W~H 
for low Higgs masses, when H bb is the dominant decay mode, and e'^e~ — t- Hbb, when 
the Higgs decays resonantly to W^W~ with a large branching fraction. The process 
e^e~ — )■ W^W~ Z with Z decaying to bb is less prominent, since the Z ^ bb branching 
fraction is only about 15%. We do not compute these additional three-particle processes 
in the present paper, since they are further suppressed when invariant-mass cuts in the 
top anti-top signal region are applied. Moreover, in a realistic analysis these contributions 
would be eliminated anyway by invariant-mass cuts on the resonance (Higgs or Z) decay 
products. 

3 Computation of the four-electron matching coeffi- 
cients at NLO 

Ik) 

The four-electron matching coefficients C\J originate from the hard contributions of the 
e^e~ forward-scattering amplitude. The hard momentum region expansion dictates that 
the top-quark self-energy insertions are treated perturbatively, since the top lines are 
formally far off-shell, p\ — ml ~ 0{mf) ^ S(p^). Accordingly the calculation of the 
coefficients C^g"* is performed in fixed-order perturbation theory in the full electroweak 
theory with no resummation of self-energy insertions in the top-quark propagator [M] 
and supplemented with an expansion of the amplitudes in 6 = s/{4mf) — 1 ~ u ^. The 
corresponding contributions to the imaginary part of the forward-scattering amplitude 
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are then reproduced by those imaginary parts of the short- distance coefficients CiJ 
which can be identified with the W^W~bb final state. The extraction of these imaginary 
parts is more conveniently performed by evaluating the relevant cuts of the hard e~^e~ 
forward-scattering diagrams directly. The leading contribution arises from the cut one- 
loop diagram with a top and anti-top propagator; however, the hard region of this 
diagram vanishes in dimensional regularization to all orders in the ^-expansion. 

A non-zero imaginary contribution to the short- distance coefficients Cjg'' arises, in the 
unitary gauge, from the cut two-loop diagrams shown in Figure [T] expanded near thresh- 
old (in 6 = s/lAm'f) — 1). Including cuts of the top lines is consistent with the expansion 
of the imaginary part of the resummed top-quark propagator in the hard region |34j^ . 
The contribution to the W^W~bb cross section from diagrams hi-hio can be interpreted 
as the bW~ pair originating from a nearly on-shell anti-top decay, while the bW^ pair is 
produced non-resonantly, either from a highly virtual top (diagrams hi-h^), or without 
an intermediate top as in the truly single-resonant diagrams h^-hiQ. The contributions 
from t bW~ cuts, which are not displayed in Figure [H yield identical results, by virtue 
of CP- invariance. The calculation of the amplitudes in Figure [1] amounts to the calcula- 
tion of the squared and phase-space integrated matrix element of the on-shell processes 
e^e~ — )■ t W~b and e'^e" — i W'^b in ordinary perturbation theory with no width added 
to top propagators and divergences regularized dimensionally for consistency with the 
calculation of the resonant contribution to ([2]). The result is parametrically of order 
a^^, and hence represents an NLO correction relative to the 0{a\-^v) leading-order 
cross section. 

Beyond the NLO non-resonant terms that we calculate in this paper, the NNLO 
contributions are given by the 3-loop cut diagrams containing the a^-corrections to the 
amplitudes of Figure [H while the next order in the 5-expansion of the two-loop diagrams 
contributes only at N^LO. It is also interesting to note that pure background diagrams, 
i.e. diagrams without any top lines, start to contribute to the imaginary parts of the 
short-distance coefficients C^^ at N^LO through the W^W~bb cuts of pure EW 3-loop 
diagrams. An exception to this are the three-particle final states with resonant decay 
of one of the particles to W^W~bb with large branching fraction as discussed in the 
previous section. 

Turning now to the calculation, it proves convenient to cast the contributions from 
diagrams hi-hiQ (at order 5° and for mi, = 0) in the form 



^Recall that applying the Cutkosky rule for an internal line leads to the replacement of the propagator 
by its imaginary part. The leading term obtained by the expansion of the imaginary part in the hard 
region leads to a Dirac delta function, which is interpreted as if a top line with no self-energy insertion 
is cut. 




(10) 



where we have introduced the variable 
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which corresponds to the invariant mass of the bW~^ system formed by the b and 
hnes that are cut. For diagrams hi-h^, where the bW~^ system comes from the decay 
of the top hne, the momentum pt corresponds to the momentum of the top hne. The 
upper hmit on the integration in ( ITOl) corresponds to the kinematic hmit associated 
to the ibW^ final state when 5—7-0. The kinematic lower limit is given by = M^. 
The case of loose cuts on the bW^ and bW~ invariant masses that are discussed in 
Section H] can be incorporated easily in the computation of the four-electron matching 
coefficients by setting the lower bound of the integration to = ml — A^, with the 
parameter A = A_ defined in ( IT9l) below. The factor [ml — p2^('^~3)/2 (fipl) regularizes 
the endpoint singularity in diagram hi at pi = due to the two top propagators that 
make Hi behave as 

M^^^l ^ const x^^ 2^. 12 

The use of dimensional regularization within the threshold expansion thus provides a 
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natural way to regularize the phase-space singularities in the hard region where the top 
propagators have no width. Dimensional regularization must be used for consistency, 
when the pure QCD calculation is performed with the dimensionally regularized non- 
relativistic effective theory. In the calculation of the other diagrams h2-hiQ we can set 
d = 4, since they contain at most one internal top propagator leading to an integrable 
square- root endpoint singularity in ( fTOl) . 

The result of the calculation can be written as 



r(l) 



327r2a2 pBom 

s irtt 



QlC^^{s) - 2vtQtC^z{s) + Czz{s)v-, 



6^2 



mt 



7r2 A 



+ C^^(s) [QlhX + 2Qi,Qthl + 2Qthl + Qlh^ + 2Q^,h^ + hio] 



- 2QtVthX + 2QtathX^ - 2{QbVt + Qt{vb + ab))h. 



2atQbh^ - 2vthl + 2ath^ - 2—Qthl 



2Qb{vb + ab)h8 - 2(vb + ab + Qb—jhg - 2—h 



no 



Czz(s) 



v'^h\ + aih^ - 2vtathX^ + 2{vb + ab){vthl + ath 



+ 2^{vthX - athf) + {vb + abfhs + 2^{vb + ab)h, + 



+ CJs) 



QthX + /i6 + Qbhj 



+ Cz{s) VthX + ath^ + —he + {vb + ab)hY 

Sin 



1 , 



(13) 



where the functions 

Tn?' 

C^^{s) = — — 



C,z{s) 



QeVefn 



4{s - Ml 



s^,, s 



Cz(s) 



Czzis) 



slis - M|) 



S jVe + CilWt 

4(s - M|)2 



retain the exact s-dependence from the photon and Z-boson propagators, and 



.Born a\Vtb\'^mt (1 - x)2(l + 2x) 



16s 



X 



(14) 



(15) 



with X = M^/mf is the tree- level top decay width obtained from t — ?■ bW'^ with the 
bottom-quark mass set to zero. The dimensionless functions 



hf ^hfix,A'/m', 



(16) 
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Figure 2: NLO non-resonant contributions to the cross sections computed at s = 4mf as a 
function of (A/m^)^, and with values for the input parameters given in fl20p . The maximal 
value of (A/m^)^ corresponding to the total cross section equals 1 — M^/m^ = 0.7815. 
The solid (black) line is the total contribution cTnon-res in (fT3|) . and the contributions 
from the diagram hi, from the sum of diagrams h2-h^ and from the truly single- resonant 
diagrams h^-hiQ are shown by the dashed (purple), dash-dotted (red) and dotted (blue) 
lines, respectively. 



with X = V,A,VA, are obtained from diagrams i = 1,...10, and depend only on the 
ratio^ {Myy/rrity and {^/mtY but not on the centre-of-mass energy y/s. To define 
the function h\ we subtracted from the integrand of diagram hi its leading singular 
behaviour as pi mf (see (l33i) in Appendix |A]). Adding back the subtracted expression 
yields the terms in the first line of ( 1T3|) upon pf integration. Integral representations of 
the functions hf which can easily be evaluated numerically are provided in the Appendix. 

At NLO we could put s = 4m^ in the Cj(s) functions appearing in (IT^ and in the 
overall 1/s factor; keeping the exact energy dependence slightly improves the quality 
of the effective field theory (EFT) expansion. Since this energy dependence is mild 
compared to that of the resonant contributions within the range of energies relevant for 
the threshold region, the non-resonant contribution aion-res results in an almost constant 
shift of the cross section. 

The size of the NLO non-resonant contributions to the cross section evaluated at 
s = Amf is shown in Figure [2] as a function of A'^/mf. In this Figure the total cross 
section corresponds to the maximal value of A^/m^ = 1 — M'^/m'l = 0.7815. The 

^Here we already allow for an invariant- mass cut on as discussed in Section |4l With no cut, 
A = Mw and A ^ y/m^ - M^. in 
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analysis reveals that the total (solid line in Figure E]) is negativ^ and follows closely the 
contribution from diagram hi alone (dashed line) for all values of A^. The sum of the 
contributions of diagrams /i2^^4 (dash-dotted line), where the bW~^ is produced from a 
highly virtual top quark, is rather small, less than 3 fb for all values of A^. The individual 
contributions from diagrams /13 and are comparable in size with that of hi, but they 
cancel to a large extent. The sum of the truly single-resonant diagrams h^-hio (dotted 
line in Figure [2]) is positive and also small, reaching at most 7 fb for the total cross 
section {i.e. when A^ = rnf — M^). The smallness of the single-resonant contributions 
is due to large cancellations among diagrams /15, h^, and hiQ which individually become 
more important than diagram hi for large values of A^. 

For easy use we provide a compact fit formula for the NLO non-resonant contributions 
to the total cross section (corresponding to A'^/mf = 1 — x) in a reasonable range of top- 
mass values. In the range mt = 160-180 GeV, a quadratic fit to the functions hf[x,x), 

hf'^\x,x) =a + b6„,,+c6'^^ , 6mt = ^^J^^y " ^ ^ i^V 

is accurate to a precision better than one per mille. The coefficients a, b, c are given 
in Tabled! When the functions hf'^^{x,x) are combined into the total non-resonant 
contribution (TnoLros according to ( |T3l) . the precision of the quadratic fit compared to 
an exact (numerical) evaluation is better than 5 ■ 10~^ for 160 GeV < rrit < 180 GeV. 
Omitting the quadratic term in hf'^^{x,x), i.e. setting the c-coefficients to zero, the 
precision of crion-rcs is still about one per mille or better. 



4 Invariant- mass cuts in the effective theory 

The measurement of the ti threshold shape under realistic experimental conditions will 
require to apply cuts in the kinematic variables of the observed particles in order to select 
W~^W~bb events. For a general treatment of such phase-space cuts the EFT framework 
will have to be modified, and the bottom-quark and ly-boson degrees of freedom and 
eventually those of their decay products have to be included. However, cuts on kinematic 
variables that are not sensitive to the angular or momentum distributions of the top 
decay products, such as cuts on the invariant masses of the bW~^ and bW~ subsystems, 
can be implemented in the unstable-particle effective-theory calculation, as discussed 
in [35ll36f37j . 

We consider symmetric cuts on the invariant masses of the top (Mt) and anti-top 
(Mf) decay products of the form 

mt - AMt < Mtf <mt + AMt , (18) 

•^The non-resonant cross section can be negative, since it is defined by analytic continuation to four 
dimensions of an expression that is divergent without regularization. Only the sum of resonant and 
non-resonant terms in ([2]) represents a physical quantity. 
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hY 
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-0 01716 
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-0 005359 




0.005727 


0.01925 


0.02413 




0.07369 


0.2965 


0.3855 


h'j 


0.001891 


0.006638 


0.004723 


hs 


-0.007086 


-0.01070 


0.009648 


kg 


-0.02321 


-0.06746 


-0.02936 


hio 


-0.5453 


-2.282 


-3.081 



Table 1: Coefficients for the ffi functions /i^ ' (x, x). 

or equivalently 

-A2_<M5--m?<A2 , (19) 

with = 2mt/S.Mt =F AM^^. Invariant-mass cuts on the top decay products translate 
to cuts on the momenta of the top (pt) and anti-top [pi) circulating in the potential and 
hard loops of the e^e~ forward-scattering amplitude. The cuts of the form (lT9l) can be 
implemented by inserting a product of step-functions 6'(A^ — [pi — 'm1])6{p1 — -|- A^), 
and similarly for pf, in the cut loop integral. In the effective theory the same step 
functions have to be introduced into the loops in the matrix elements of the resonant 
contributions, which reproduce the potential region of the full-theory diagrams. In the 
hard region, these step functions modify the calculation of the matching coefficients of 
the four-electron operators, which then become dependent on the invariant-mass cuts 
A^. At NLO this dependence enters explicitly through the lower limit in the integration 
that defines the hf functions, see ( ITO|) . (fT6|) . At this order, no dependence on the upper 
invariant-mass cut A^ appears because the phase-space integration of the ibW~^ final 
state in the hard diagrams of Figure [T] sets an upper kinematic limit p^ < ; hence we 
define A = A„. 

The two limiting cases for the invariant- mass cut corresponding to a loose cut, 
A^ ~ mf (corresponding to AMt ~ mj), and a tight cut, A^ ~ mtTt (corresponding to 
AMf ~ Ft), yield particularly simple prescriptions for implementing invariant-mass cuts 
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in unstable-particle effective theory, as it has been discussed in the context of ly-pair 
production near threshold [35]. For tight cuts the matching coefficients of the four- 
electron operators vanish because in the hard region we have — mf\ ^ A^. Then 
the expansion of the product of step functions in small variables provides the factor 
^{~\Pt ~ ''^t I) — 0- Therefore four-electron operators do not contribute to the cross sec- 
tion for tight invariant-mass cuts. The same does not apply for the contributions from 
the resonant part, because there — = 2mtrQ — f^r^ A^, since the residual momen- 
tum r = (ro,r) has the characteristic scaling of the potential region r ~ {rritv'^ , rritv) . 
The situation is reversed if we consider loose cuts. In this case the implementation of 
invariant-mass cuts on the loop integrals from the effective theory has no effect at the 
leading order because we can drop the momentum dependence from the step function 
following the assumption A^ S> mtro, r ^ ~ {rritvY, and we are left with 6'(A^)6'(A^) = 1. 
In particular, this means that the NLO QCD interference contributions are not affected 
by loose cuts and thus vanish, as it happens for the total cross section. The latter is 
confirmed by the calculation of the QCD interference effect with an invariant mass cut 
satisfying A^ ^ mtTt [27] , which yields a correction of order m^Fj / relative to the LO 
resonant term, which is indeed parametrically suppressed relative to a NLO correction 
of order a/oew when A^ ~ m^. However, as already mentioned, the loose cut needs to 
be taken into account in the calculation of the matching coefficients of the four-electron 
operators as done in Section [31 

It is interesting to compare the above procedure for the treatment of invariant-mass 
cuts in the unstable-particle effective theory with the phase-space matching approach 
introduced in the recent paper [37]. In the latter, the matching coefficients of the four- 
electron operators (called phase-space matching conditions) are determined in the follow- 
ing way. First, cuts of the form f|T8|) are introduced on the top and anti-top momenta in 
the NRQCD loop diagrams, the latter corresponding to the resonant contributions in ([2]). 
Then the resulting cut loop integrals are expanded assuming that A^ ^ (mtu)^,mtFj, 
which is formally correct as long as A^ ^ m^, since NRQCD describes the non-relativistic 
tt configurations and A^ acts effectively as a cut on the non-relativistic momentum 
r*^. The A-dependent terms that result from this expansion yield a series of the form 

Ff/A X X]nm[("^<-'^i/^^)" ^ (^^/■'^t)"']! with n, m = 0, 1 Powers oi A?/ml enter first 

from cut diagrams with NNLO relativistic corrections that come with a factor f'^/mf, 
and are a signal of the breakdown of the series for A ~ m-f. The whole computation is 
equivalent to the non-relativistic expansion of the full-theory squared matrix elements 
containing the double-resonant diagrams for e~^e~ ti W~^W~bb and their interfer- 
ence with the diagrams for e~^e~ — W~^W~bb having only either the top or the anti-top 
in intermediate stages. Contrary to the method presented in the present paper, the 
non-resonant contributions to the full-theory matrix element coming from the square of 
single-top and pure background diagrams (the so-called remainder contributions in [37] ) 
cannot be determined within the phase-space matching approach. The path taken by 
the authors of [37] is to compute the remainder contributions by comparing the NRQCD 
cross section with the relativistic Born cross section for e~^e~ — )■ W~^W~bb as obtained 
with MadGraph/MadEvent [H]. The comparison yielded that the remainder terms are 
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below 5 fb for invariant-mass cuts AMj < 35 GeV, and hence were neglected given the 
theoretical precision goal aimed at in |37]. This approach is currently not feasible at 
NNLO, where no external tools exist to determine the remainder terms, and the a^- 
correction to the coefficient functions of the four-electron operators must be determined 
by the methods described in the present paper. Also, the present method allows us to 
consider very loose cuts with A ~ up to the total cross section. On the other hand, 
assuming that there is no kinematic or dynamical enhancement in the QCD corrections 
to the remainder contributions, part of the a^-corrections to the non-resonant contribu- 
tions has already been analyzed in [37], again for intermediate-sized cuts, which in our 
approach correspond to NNLO contributions not yet known. 

We have checked that the expansion of our result (fT3l) for the non-resonant contri- 
butions in A? /ml <C 1 matches the series obtained in [37] when = for the Fj/A and 
TtA/m^ terms. In this limit, diagrams h^-hiQ, which are part of the remainder contri- 
butions not calculable within the phase-space matching approach, first contribute with 
TtA^/mf terms. 

4.1 Comparison to the e+e~ W^W~bb Born cross section 
from Mad Graph 

Before discussing the final result for the W~^W~bb cross section in the top anti-top reso- 
nance region including Coulomb resummation, we compare the EFT prediction without 
strong-interaction effects to the full-theory (Standard Model) W~^W~bb Born cross sec- 
tion computed with MadGraph/MadEvent [H]. The input parameters are chosen to 
be 

Mz = 91.1876 GeV , Mw = 80.398 GeV , Mh = 120 GeV , 

= 1.166367 X 10"^ GeV'^ , mj = 172.0 GeV , Va = I , (20) 

whereas the on-shell Weinberg angle = M-^/Mz and the fine-structure constant in the 
G^-scheme, a = v^G^A-f^s^/vr, are derived quantities. Tree-level amplitudes generated 
by MadGraph implement the fixed-width prescription for the top-quark propagator. For 
the numerical value of Fj we use the Born formula (IT^ in both MadGraph and the EFT 
predictions, thus neglecting a^-corrections to the top-width value. MadGraph uses a 
finite bottom-quark mass whereas we neglect it in the EFT calculation. 

In Figure [3] we analyze the invariant-mass cut dependence of the e^e~ — )■ W^W~bb 
Born cross section calculated at s = 4m^. The values shown there range from the AMj 
value corresponding to A^ = mjF^ up to the maximum value allowed by the kinematics, 
AMt^max = fTit — Mw- The (red) triangle-shaped dots represent the Standard-Model 
e^e~ — > W^W~bb Born cross section computed with MadGraph, while for the (blue) 
circular dots we have omitted in the MadGraph amplitude generation the diagrams with 
internal Higgs lines. The shaded bands represent the statistical uncertainties on the 
MadEvent integration, calculated as the cross section divided by the square root of the 
number of selected events that pass the invariant-mass cut. The small difference between 
both sets of points for large values of AM^ is explained below. The lower solid (blue) 
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A Mt (GeV) 

Figure 3: Invariant-mass cut dependence of the e^e" — >■ W^W~bb Born cross section at 
s = Am"^. The AMt values range from AA/^ — 0.734 GeV corresponding to — TrifTf to 
the kinematic bound, AM^ max = rrit — Mw = 91.602 GeV. The (red) triangle-shaped dots 
represent the Standard-Model cross section computed with MadGraph; for the (blue) 
circular dots we have omitted the diagrams with internal Higgs lines. The corresponding 
shaded bands show the statistical uncertainty of the MadEvent integration. Lower solid 
(blue) line: NNLO EFT approximation given in f l2T|) . Upper solid (brown) line: LO EFT 
tight-cut prediction f pOj) . Dashed (brown) line: = NNLO results with invariant-mass 
cuts from [37] (not including the remainder contributions). The LO Born EFT result, 
which does not depend on AMt, is not shown in the Figure. It equals cr^^^ = 0.1139 pb 
at s = 4mf. 



line is the NNLO EFT approximation to the full Born result, given by 

_ ^(0) , _.(!) 1^(2) (2) I ^(2) 1^(2) /9-|N 

O^Born - %,„^=o + ^non-rcs + ^t,2 + ^^P-wavc + ^bil + ^abs • {^^) 

The first term on the right-hand side is the LO cross section ([6]), with the QCD radiative 
corrections switched off {i.e. setting = 0). At s = Amf it amounts to c^^^ = 
0.1139 pb. Loose cuts on the bW~^ and bW~ invariant masses contribute first at the NLO 
in the EFT power counting through the non-resonant contributions (Jnon-res discussed in 
Section |3l The w^-suppressed S'-wave current describing the top anti-top pair production 

(2) 

in the effective theory gives the NNLO correction • Using the non-relativistic equation 
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of motion for the top quark it can be taken into account by replacing Gq\o,0;S) — )■ 
(1 - £/{3mt))G''c\0,0;£) in The term 

^p-wave is generated by the matrix element 
of the production (decay) P-wave operators, which are f-suppressed with respect the 
leading 5- wave operators It reads 



(2) 



P-wavc 



\'^p,P-wa.vc) ^ V-'p,P-w!we) 



4iVe 

3mf 



ImG 



(2) 



P-wave ' 



(22) 



with 



p, P-wave 
(a) 



47ra 



s-Ml 



^^„P-wave=-4™^3^ 



(23) 



and 



G 



(2) 

P-wavc 



Air 



^.3/2 

nit) 



(24) 



the P-wave component of the free Green function at zero distance that can be obtained 
by setting = in the MS-renormalized expression for the P-wave Coulomb Green 
function [T3] . The term a^^ originates from the effective-theory matrix elements with a 
single insertion of the NNLO corrections to the bilinear quark terms in the non-relativistic 
Lagrangian, 6C = l/(8m^) X'ilj\{d'^ + imtTtY\lJt, and the corresponding anti-quark terms. 
Written as corrections to the zero-distance Green function, it reads 



(2) 
bil 



c-p -I- C-p 



2Arjm {5G^^^ + 5Gf] + 5G 



(25) 



where 5G^ 



(2) 



and SG'^^} correspond to the insertion of the (^"^/(Sm^), iV td"^ / {Am^) 



and — r^/(8mt) operators, respectively, and read 



5G 



(2) 
kin 



(2) 



5G 



5G^'} 



32^ 



3m^ 



nit) 



Tt 



IGtt mt 
1 



mt 



1/2 



327r (-£/mt)V2 



(26) 

(27) 
(28) 



The terms ^Gp/ and SG^i are pure electroweak contributions. While the former corre- 



sponds to the lifetime-dilatation correction studied in [21], 5Gr? has not been consid- 
ered so far among the NNLO electroweak corrections. Numerically, the latter term is 
rather small. Note that can be obtained by replacing rrit — >■ nit — iTt/2 in © 

(2) 

and expanding to first order in the correction; similarly, 6G„2 follows from replacing 
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Ft —> Tt + iTf/{4mt) in ([6]) and the first-order expansion in the correction. In both 
cases, only the term with = is needed for comparison with the Born cross section. 
The last term in (12T|) . crabs' arises from the absorptive parts in the matching coefficients 



of the production operators O^^ and O^""^ that were determined in [31]. It reads 

^Ibi = 2 K'^C^"^'"''' + C^^^C^'^)'"^^] 2N, Re G^^'> , (29) 

with Gq*^^ the free Green function. The explicit expressions for the coefficients Cp^'''^^^ 
and Cp"''^"'^^^ can be found in [31] (there denoted C^'^^^ and C'^'^^^, respectively) The 
term a^j^g reproduces the interference of single-resonant and double-resonant e~^e~ — )■ 
W^W~hh diagrams in the kinematic region where both hW^ and hW~ have invariant 
masses of order m^. Within the threshold expansion [12], these contributions arise from 
precisely the cut diagrams hi-h^ in Figure [U but with one hard and one potential loop, 
the latter corresponding to the loop containing the top and anti-top propagators. Let us 
recall that for the non-resonant contributions the two loop momenta have to be taken as 
hard, so the absorptive parts in C^"*'^^^ and Cp""^'"^^^ and the non-resonant contributions 
come from different loop-momentum regions and there is no double-counting. The ad- 
ditional ^-suppression of the hard-potential with respect to the hard-hard contributions 
from diagrams hi-h^ is explained by the different power counting of the loop momentum 
in the hard (rg ~ r ~ rrit) and potential (rg ~ f'^/rrit ~ rrttv'^) regions. The sum of all 
the NNLO corrections in fl2T|) gives a positive shift of 4.1 fb for s = Arnf. 

Already for values of the invariant-mass cut AMt > 5 GeV ~ 3.5 we find an 
excellent agreement between the full-theory Born cross section computed with MadGraph 
and the effective-theory prediction for loose cuts. To account for invariant-mass cuts in 
the range AMf < F^ we need to implement the tight-cut prescription in the effective- 
theory calculation. This is done in the upper solid curve of Figure [3l which is obtained by 
introducing the proper step functions into the one-loop cut integral of the LO effective- 
theory result. Precisely, the LO tight-cut prediction reads 



^tt, tight 



Op -I- c-p 



2iVJmGi°\g,,, (30) 



with 



j^^(o) _ r? f d^r 9{AM, + f + ro - £:) ^(AM, + f - ro - £) 



'0, tight 



(2vr) 



- 2 ^ ' 2mt) ^4 



f-'-o-£)= + i 



. !^ (, , iL)"V ^ ,.,ta„ { ^) . h.c.) . (31) 

which we evaluate a.t E = for the result displayed in Figure |3i We see that this 
expression agrees well with the full Born cross section for very small AMj. The agreement 

^The imaginary part of the field-rescaling factor called w in |33|34j is taken into account in the 
calculation of these coefficients. 
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could be extended to larger values by implementing the cuts in the NNLO terms, but 
we do not pursue this here. 

Finally, the dashed curve corresponds to the = NNLO formula with invariant- 
mass cuts obtained in [37] (not including the remainder contributions), which contains 
the same NNLO pieces as our result (1211) but for 6Gj-2. As mentioned before, the pieces 
depending on the invariant-mass cut in the formula in [37] match the leading and sub- 
leading terms of the expansion in A^/mf of our NLO non- resonant contributions, and the 
two results should therefore agree well when AMj ^ AMj max- Figure |3] demonstrates 
this and also shows that as expected, the difference between our result (solid blue line) 
and the formula from [37] increases with AMf. at AMt = 35 GeV it amounts to about 
4 fb, while for the total cross section (AMf max) it reaches almost 9 fb, which roughly 
agrees with the size of the single-top non-resonant contributions, that are not accounted 
for in the phase-space matching approach. 

As mentioned above. Figure [3] shows two sets of MadGraph points. The red triangles 
represent the complete tree-level cross section. They agree with our effective-theory 
prediction for loose cuts AMj > 5 GeV, but for AMf > 40 GeV and for the total 
cross section they rise to values of up to 3 fb above our result. This small, but visible 
difference is mostly due to diagrams with intermediate Higgs lines which account for 2 fb 
of this difference. Eliminating the Higgs diagrams from the MadGraph amplitude gives 
the blue dots which agree with the effective-theory result also for large AM^ within the 
statistical uncertainty of the MadEvent integration. The Higgs contribution originates 
almost exclusively from phase-space regions where the bb invariant mass is close to the 
Higgs-mass value Mh = 120 GeV, i.e. from on-shell Higgs bosons decaying into a bb 
pair. Despite the smallness of the bottom Yukawa coupling, such diagrams can yield 
sizable contributions because the Higgs width is small and the branching fraction of 
H ^ bb is large (about 70%). The distribution of the Higgs contributions with respect 
to the invariant masses of the bW~^ or bW~ pairs is relatively broad, while the same 
invariant-mass distributions of the contributions with (anti)top lines are sharply peaked 
around the top mass. This explains why the Higgs contributions are only visible for very 
loose or no cuts. Such small contributions from Higgs diagrams constitute a reducible 
background which can easily be eliminated by cuts on the bb invariant mass when one 
is only interested in the tt resonance region. This justifies that we have not considered 
contributions from nearly on-shell Higgs bosons in our effective-theory approach although 
they formally count as NLO contributions as explained in Section 12.21 

The energy dependence of the effective-theory prediction fl2T]) for the total cross sec- 
tion (upper panel), AMt = 15 GeV (middle panel) and AMt = 5 GeV (lower panel) 
is shown in Figure S] in the interval -y/i = (334, 354) GeV. The blue dots correspond 
to the full-theory Born W~^W~bb cross section, while the solid (blue) line is the corre- 
sponding prediction at NNLO from the unstable-particle effective-theory. For the total 
cross section the tiny difference between both results is mainly due to the Higgs contri- 
butions explained above. The agreement improves for smaller values of AMt, as seen in 
the middle panel of Figure |H For AMt = 5 GeV the effective-theory expansion quickly 
degrades as we move away from the threshold. This is evident below threshold, where 
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Figure 4: NNLO EFT approximation to the Born W~^W~bb cross section with (sohd 
blue hne) and without (dotted black line) the NLO non-resonant contributions as a 
function of the centre-of-mass energy. Upper, middle and lower panels correspond to 
AMt = (AMf max, 15, 5) GeV. The (blue) dots are the full-theory result computed with 
MadGraph, and the shaded band is the statistical uncertainty of the MadEvent integra- 
tion. 
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the effective-theory result becomes negative. The large negative shift of about 24-39 fb 
for the total cross section given by the NLO non-resonant contributions is clearly seen 
in Figure H] once we draw the NNLO effective-theory result without the non-resonant 
contributions (dotted line). The shift becomes larger as we tighten the invariant-mass 
cut, up to 38-48 fb for AMt = 15 GeV. 

5 Final results 

In this section we compare the size of the NLO non-resonant and QED corrections to the 
LO effective-theory approximation ([6]) for the e^e~ — )■ W^W~bb cross section. This now 
includes the summation of Coulomb corrections proportional to (a^/f)" to all orders in 
the strong coupling, which produces the characteristic peak structure in the top anti-top 
resonance region. As mentioned in Section 12. 2[ we do not discuss the pure NLO QCD 
corrections, which have already been studied in the literature. Moreover it is known 
that the NLO QCD corrections are rather large (and negative) and cancelled to a large 
extent by the NNLO QCD corrections p[7|H|9pi]|llfl2] . Hence the conclusions that can 
be drawn by comparing the electroweak corrections with the NLO QCD corrections may 
not be valid once higher-order QCD corrections are considered. For the results in this 
section we choose 0^(30 GeV) = 0.142 for the value of the QCD coupling that enters in 
the LO Coulomb Green function. The relative sizes of the NLO electroweak corrections 
with respect to the LO result are displayed in Figure [5l 

The upper (blue) solid line shows the ratio o'qej^/o'I^'' (in percent), where the QED 

correction o"nln is obtained as 



The QED contribution represents a correction of about 2% above threshold and rises 
to a maximum of 7% just below the peak. The lower (red) solid line is obtained as 
the ratio crion-rcs/cr^^f ^ (in percent) with AMt set by the kinematic bound. The non- 
resonant contributions give a constant negative shift of about 3% above threshold. Below 
threshold the relative size of the non-resonant corrections is very large, since the LO result 
rapidly vanishes, reaching up to 19%. Hence below threshold they represent the leading 
electroweak correction to the total ti cross section. The sum of both NLO electroweak 
corrections compared to the LO cross section is shown in the middle solid (black) line. 
We observe a partial cancellation of the QED and non-resonant corrections in the peak 
region and at energies above. A sensitivity to the invariant-mass cut AMt in the bW~^ and 
bW~ subsystem enters first at NLO through the non-resonant contributions. Restricting 
the available phase-space for the final-state particles by tightening the invariant-mass 
cuts AMt makes the non-resonant contributions even more important. This is shown by 
the dashed lines in Figure corresponding to AMt = 15 GeV, for the relative size of 
the non-resonant correction (lower dashed line) and of the sum of the two electroweak 
corrections (middle dashed line). 
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Figure 5: Relative sizes of the QED and non-resonant corrections with respect to the tt 
LO cross section in percent: o'QED/^tf'' (upper sohd blue line), anJn-ves / cr^'^^ for the total 
cross section (lower solid red line) and AMj = 15 GeV (lower dashed red line). The 
relative size of the sum of the QED and non-resonant corrections is represented by the 
middle (black) lines, for AMj max (solid) and AMf = 15 GeV (dashed). 



Aside from the pure QCD corrections, the effective-theory NLO prediction for the 
e'^e~ — )■ W~^W~bb cross section is displayed by the solid lines in Figure [61 The upper 
panel corresponds to the total cross section and for the lower panel AMf = 15 GeV. The 
absolute size of the non-resonant correction is given by the difference between the dashed 
lines, which only include the QED NLO correction, and the solid ones. These negative 
shifts amount to 27-35 fb for the total cross section and 39-46 fb for AMt = 15 GeV for 
a/s in the interval (338, 350) GeV. 

6 Summary 

We presented an analysis of the e~^e~ — W^W^hh process in the top anti-top resonance 
production region ^/s ^ 2mt extending the methods developed in [Mf35] for VT-pair 
production to the more complicated case of top quarks. Our result may be considered as 
the first complete NLO calculation of top-quark pair production near threshold, when the 
top width is accounted for, which necessarily requires to consider the final state W^W~bb 
for consistency. The result is valid for the total cross section and with invariant-mass 
cuts on the Wb systems as long as the cut AMt is significantly larger than the top width. 
We also included the case of tight cuts in the comparison with the full-theory Born cross 
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Figure 6: e'^e — ?■ W^W bb cross section with LO QCD effects and NLO electroweak 
corrections at energies close to tliresliold. Tlie daslied (blue) lines represent a^^^ + cTq^^ 

and the solid lines a[^^ + Cqeq + crion-rcs- The results in the upper panel correspond to 
the total cross section (AMt = AMt^max), while for the lower panel AMj = 15 GeV. 



section. 

We find that the non-resonant contributions reduce the cross section in the top anti- 
top resonance peak region by up to 5%, which is relevant in view of the high-precision 
top mass and width determination anticipated at the ILC. In the energy region below the 
peak the non-resonant correction is significantly larger, since the resonant contribution 
rapidly decreases with energy, while the non-resonant one is nearly energy-independent. 
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These conclusions agree with the recent work |37j. The approach pursued in the present 
paper allows us to extend them to the case of the total cross section and rather loose 
cuts, and it also determines the full non-resonant background within the effective theory 
approach rather than using external information. The final numerical result turns out 
to be dominated by the off-shell contribution from the double-resonant diagram hi, but 
only due to a large cancellation within classes of other diagrams. 

Our work provides a step towards predictions including consistently the top width 
in effective theory calculations of ti production, more precisely W~^W~bb production, 
near threshold also in higher orders in the expansion in a^, ctEw and 5 ~ f^. Going 
to higher orders will eventually be required to cancel the finite-width divergences that 
appear in the NNLO and N'^LO pure QCD calculations that are already complete or 
near completion. 
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A Formulae for the hf functions 

We present here integral representations of the functions hf{x,y) defined in f|T6|) and 
used in f|T3l) . Their arguments are x = / and the cut-dependent integration limit 
y = A^/m^ = 1 — A'^/mf satisfying x < y < 1. 



Diagram hi. 

hi{x,y) - 



with the "off-shell" top width defined as 

, 2 



TAt) = rf°'"^A- ^^-^ ^ . (34) 

The last term within the square brackets in the first line of fl33l) is a subtraction term 
which regularizes the integrand at t = 1. The integration of minus this subtraction 
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term using dimensional regularization yields the terms quoted in the first line on the 
right-hand side of f lT3|) . 

Diagram h2'- 



1 



hY(x,y) = -— -7-- r / dt 

2^ '^^ 47r2(a;-l)2(2x + l) 4 

t (2x2 + 7x - 5) - X (2x2 + 15x + 31) , , , 
^—j. /+(^' ^) 



+ y Y — ^ 0^ + + 4) - (2x2 + 12x + 53) - tx(24x + 85) - Gx^) 

(35) 



1 



hr, (x, y) = — — -—, / dt 

2^ '^^ 47r2(x - l)2(2x + 1) ' 



y 



t (2x2 + 5x - 3) - X (2x2 + 13x + 17) , , , 
Y—j. /+(^' ^) 



+ y ^ — ^ (^t^ + t3(x + 2) - ^2 (2x2 + 14x + 27) - tx(20x + 59) + 6x2^ 

(36) 



with 

. m f + . (37) 

I v/(9 - t)(l - t)(t - x) - t(t - 5) ± x(t + 3) ' ^ ^ 



Diagram /13: 

^ = - 

^7r2(x - l)2x(2x + 1) 



1 

hl{x,y) = —- —r— — — / dt 



I (t^x + ^2 (4x2 - 6x + 2) + t (-2x^ - 29x2 + 2) 
+ X (4x^ + 36x2 ^ 3, _ ^3^\ jj^^ ^) 



+ rLJ. ( + t3 _ 4^ _ 4) + ^2 /_23,3 _ 283,2 ^ ^ 4N 

8t V 1 — ^ ^ 

- tx (24x2 + 69x + 92) - 6x=') I , (38) 
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(t^x + (4x^ -6x + 2) +t (-2x^ - 23x^ + 20x + 2) 
+ X (4x^ + 24a;2 - 33x + 7) ) /„ {t, x) 



^— ^ W ^ — - (t^x + t^{x - 6)x + (-2x3 - SOx^ + 53x + 16) 
+ tx (-20x2 - 75x + 56) + 6x' 



Diagram /i4: 



^^"^^"'^^ = - 64vr^(l-!:)^(2x + l) / ^^T^ 
t(2tx + t- 17x - 1) + 23x 



(t-x) 



6tx^ 



t^x + (2x2 - 5x + 2) + t(x + 2) - x(86x + 13) , , , 

H , ^= /_ (t, x) 

3v/(9-t)(l -t) x2 

.2 N - 4^3 - 106^2 + 604t - 495 
+ - (t^ - 5t2 - 25t + 93) g^^>(t, x + , g^^>(t, 

+ - {-t^ + 21^2 - 119t + 99) g^'\t,x) - -((9 - t)(l - t)) (^^^^^^a;) \ 



3x /" 9 — t 

^^^^^y^ = 64.^{i-xn2x+i)l ^'VT^ 

t^x - 2t(4x + 1) - 7x 



^ ' 3tx2 

+ (2x2 - 5x + 2) + t(21x + 2) + (7 - 30x)x , , , 



3^/(9 -t)(l-t) x2 

1/ X/ x2 mw X -t^ + 4t3 + 74t2 - 284t + 207 .ix, , 
- t - 3 t - 1 2^(0) t,x + , 9^ Ht,x) 
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with 



and 



+ 1 (t3 _ 2it2 + ii9t _ 99) gi^) (t, a;) + 1 (^(9 _ t) (1 _ t)) ^(3) (t, ^) I , (41) 
9^'\t, x) = - dz , , , , In , , , , , : (42) 



2 j_i \/az^'^^l)z'^^ \d — ez + kV az'^ — bz + 



a = 4(9-t)(l -x), 

b = A^/{9-t){l-t) (x - 2)(t + 2x - 3) , 

c = t2 (2,2 + 4) _ 2t (Sx^ - 8x + 12) + 25x2 - 48x + 36 , 

d = -2t^ + t((x - l)x + 6) - x(x + 3) , 

e = -\/(9-t)(l -t) (t(x - 2) + x) , 

k = t-x. (43) 



Diagram /is: 



= 1024.^(1 -x)'(2x + l) / " * V(9 



^ ^ ^ ^ ft^x(8x + 1) + (Tx^ + 39x + 4) 

6t2a;3(t2a; + t(4_6x)+x(4x-3)) V ^ ' \ ) 

+ (-32x3 + 53x2 ^ ^ g4) ^ (283,2 ^ _ 27) + 3a;2(43, _ 3) j 

4((t-9)x2 + (t-ll)x-6) ^ ^ ^ 4(t + 3) .01. ^ 
3xV(9-t)(l-t) ^ y V , ; 



_ 4V(9 t)(l _ 4^ (^2 _ 8^ ^ ^5) ^.(0)^^^^) 

X 

- 8xit - 4)y/{9-t){l-t)j'''\t,x) - 4x(9 - t)(l - t) j(2)(t^a;) 
+ 8 ^'^ + ^(2-4x)-x-6 ^(0) ^) ^ 8 ^(9_i)(i_i) ^ (^^ ^) 

X X 



(44) 



with 



.(i), \-fT_^\ rl i (C3C^1 + Cl4)(c2 + 8C1C3) - 6CiC2C3C?2 

3 {t,x)-yi 2ciC3(ci-4ciC3)2 
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dz * ^'^'^ ~^ 2ciC3)rf2 - 3c2(cirf3 + csrfi) ( ^ ~ ^ k^/az'^ -hz + c\ 
7-1 — AciCzY \/ az"^ — bz + c \d — ez + k\/ az"^ — hz + cj 

(45) 



j»(t,x) = -(l-|) / dzz 



t' Ci(c|-4CiC3) 



' ^^^i 203^2 - C2(i3 1 f d-ez- ky/az^ -bz + ^ 

_i cl-AciCs .Jaz"^ -bz + c \d - ez + k\/ az^ - bz + c 



and 



x((t-6)t + 4x -3) 
c, = ^ 1 + 



C2 



C3 



45'(t(x - 2) + x) - t2(x - 4) - tx{x + 3) + 



d 



3 



x(2?-t)(25'-x) 

{t - 7){t - 3)^2 + (t(t + 2) + 9)x^ + {t{t + 2) - 27)te 

42'((t - 3)t2 + (t + 3)x2 + (t - 3)te) - t(t2(3x + 4) + te(3x - 23) + 13x2) 

A'P (i^ + tx + x^) - 6ztx{t + x) + ?>t^x^ 



7 t(t + X - 5) + 3x 7 zit + x) 

c^2 = — ^ , 4 = X ^ . (47) 

Diagram h^: 

1 

h^i^^y) = -TT^i n^^T^T— TT / ^^V(9-t)(l-t) 

1287r^(x — lYx^[2x + 1) y,y 

— -^^^^ ^ ^ ( t'x (6x2 + 7^ + 1) 

{t^x + t(4 - 6x) + x(4x - 3)) V ^ ^ 

+ (-21x^ + 9x2 ^ + 4) + ^2 (-40x^ - 224x^ + 36x2 _ + 84) 
+ tx (60x^ - 33x2 ^ 25x - 27) + 3x2 (4x _ g^,^ 

27 



+ ^ f- t, x) 

V(9-t)(l-t) 

- 6(t + 3)x^{3x + 2) x) + 6^/{9-t){l-t) x'^{{t + 6)x + 2) g^^\t, x) 
-6{9-t){l-t)x^g^'^\t,x) 

- 6x^ {t^x{x + 4) - 2t (5x^ + 8x - 4) + 33x^ + 20x - 24) 7°) (t, x) 
+ 12^/(9 -t){l-t) x'^{-2{t - 2)x + Tx^ - 4) j^^^^^a;) 

-6(9-t)(l -t)x^j(2)^^^^)| _ (4g) 
Diagram h-j: 



647r2(l -x)2(2x + 1) 5-t 
(t-5)(t-x)(2tx + t-15x) 

t^x - t2 (2x2 + 9x - 2) + 3t (8x2 + 3x - 2) + X (16x2 - 18x - 17) 

H , = /- (t, x) 

3v/(9-t)(l-t) x2 

4(^^x2 + t (-8x2 + X - 1) + (-4x2 + 2x + 13) x) 
3v/(9-t)(l-1y^ 

- ^ (-t^ + t2(2x + 13) - 5t(4x + 3) + 66x - 45) g^^\t, x) 
+ ^\/(9 -t)(l -t) ((t - 2)t + 28x - 39) ^(^)(t, x) 

- \{9 - t)(l - t)(t + 2x - 11) ^?(2)(t, 3;) - i((9 - t)(l - t))'/' ^7(3)(t, x) 

- ^ (t^ - t2(2x + 7) + t(20x + 3) - 66x - 45) 5^(°)(t, x) 
+ i^(9-t)(l-t) (t(t + 4) - 7(4x + 3)) fi^(i)(t, x) 

_ 1(9 _ _ _ 2x + 1) 5^(2) (t,^) + 1 ((9 _ t)(i _ t)) 3/2 ^(3)^^^^) I ^ ^4g) 
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with 



2 

and 



1 



-az^ — bz + c \d — ez + k\/ —az^ — bz + 



a = -^(9 - t){l - t) [f -2t- 16x + 1) , 

b = ^V(9-t)(l -t) ((t-6)t-8x-3)(t-2x-l), 

_ t^- 4t3(x + 3) + (4x2 + 44x + 30) - 4t (lOx^ + 35x - 9) + lOOx^ + 228x + 9 
'= 4 ' 

d = ^ {-t^ + 3t\x + 2)+t {-2x^ - 16x + 3) + x{2x + 21)) , 

e = ^\/(9 - t)(l - t) {t{-t + X + 1) + 7x) . (51) 
Diagram h^: 



(^t^(l - x) - t^(x - 5) + (3x^ + 49x2 + 44x - 72) 



(t-x)v/(9-t)(l-t) 
16t2(t(l -x) +x(x + 3)) 

- tx (2x^ + 23x2 + 107x + 144) + Gx^ (x^ + 7x + 12) ) | . (52) 

Diagram h^: 

^^^"'^^ = 8vr^(x-l)M2x + l) / 

2(t (2x3 + 7x2 - 7x + 2) - X (2x3 + 23x2 + 25x + 6)) 

(t^x - 2^2 (2x2 + 5x - 1) + t (-2x3 + 19x2 + 12x - 6) 
- X (4x3 + 16x2 - 25x + 17)) /_(t,x) 
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+ 6x2(x + 4))|. (53) 

Diagram hiQ-. 



I - (-t^ + 7t - 3) + (4t + ll)x2 + - 12x - 12^ f_{t, x) 

+ (^-^M9- t)(l-t) ^ / ^5^3 .3 47^2 + 4) 

16t2 (t2a; + t(4 _ 6a;) + xiAx - 3)) V ^ ^ 

+ {-2x^ + 19x^ - 221x^ + 64x2 + 124x + 16) 
+ (86x^ + 259x^ - 683x^ + 12x2 - 212x + 336) 
+ tx (-8x^ + 54x'^ + 201x^ - 208x2 + 148x - 108) 

+ 6x2 ^4^4 ^ ^3^3 _ ^2a;2 + _ A I (54) 



References 

[1] Tevatron Electroweak Working Group, for the CDF and DO Collaborations, 
larXiv:0903.2503i [hep-ex]. 

[2] M. Martinez and R. Miquel, Eur. Phys. J. C27, 49 (2003), |hep-ph/02073l5] . 

[3] V. S. Fadin and V. A. Khoze, JETP Lett. 46, 525 (1987). 

[4] V. S. Fadin and V. A. Khoze, Sov. J. Nucl. Phys. 48, 309 (1988). 

[5] M. J. Strassler and M. E. Peskin, Phys. Rev. D43, 1500 (1991). 

[6] A. H. Hoang and T. Teubner, Phys. Rev. D58, 114023 (1998), |hep-ph/9801397| . 

[7] K. Melnikov and A. Yelkhovsky, Nucl. Phys. B528, 59 (1998), [ hep-ph/ 9802'379] . 

[8] M. Beneke, A. Signer and V. A. Smirnov, Phys. Lett. B454, 137 (1999), 
|hep-ph/9903260l . 

[9] A. H. Hoang and T. Teubner, Phys. Rev. D60, 114027 (1999), |hep-ph/9904"468] . 

30 



[10] O. I. Yakovlev, Phys. Lett. B457, 170 (1999), |hep-ph/9808463] . 

[11] T. Nagano, A Ota and Y. Sumino, Phys. Rev. D60, 114014 (1999), 
|hep-ph/9903498| . 

[12] A. A. Penin and A. A. Pivovarov, Phys. Atom. Nud. 64, 275 (2001), 
|hep-ph/9904278| . 

[13] A. H. Hoang, A. V. Manohar, I. W. Stewart and T. Teubner, Phys. Rev. Lett. 86, 
1951 (2001), |hep-ph/0011254| . 

[14] A. H. Hoang, A. V. Manohar, L W. Stewart and T. Teubner, Phys. Rev. D65, 
014014 (2002), |h ep-ph/0107T44] . 

[15] A. H. Hoang, Phys. Rev. D 69, 034009 (2004), |hep-ph/0307376] . 

[16] A. Pineda and A. Signer, Nucl. Phys. B762, 67 (2007), |hep-ph/0607239] . 

[17] M. Beneke, Y Kiyo and K. SchuUer, Nud. Phys. B714, 67 (2005), 
|hep-ph/0501289| . 

[18] M. Beneke, Y. Kiyo and K. SchuUer, PoS RADCOR2007, 051 (2007), 
larXiv:0801. 34641 [hep-ph]. 

[19] M. Beneke and Y. Kiyo, Phys. Lett. B 668, 143 (2008), larXiv:0804.4004l [hep-ph]. 

[20] P. Marquard, J. H. Pidum, D. Seidel and M. Steinhauser, Nud. Phys. B758, 144 
(2006), [hep-ph/0607168] . 

[21] A. V. Smirnov, V. A. Smirnov and M. Stdnhauser, Phys. Lett. B668, 293 (2008), 
larXiv:0809.1927l [hep-ph]. 

[22] C. Anzai, Y. Kiyo and Y. Sumino, larXiv:0911.4"335] [hep-ph]. 

[23] A. V. Smirnov, V. A. Smirnov and M. Steinhauser, I arXiv:091 1.47421 [hep-ph]. 

[24] N. Brambi lla, A. Pineda, J. Soto and A. Vairo, Phys. Rev. D60, 091502 (1999), 
thep-ph/9903355l . 

[25] B. A. Kniehl, A. A. Penin, V. A. Smirnov and M. Steinhauser, Nucl. Phys. B635, 
357 (2002), |hep-ph/0203T66| . 

[26] A. A. Penin, V. A. Smirnov and M. Steinhauser, Nucl Phys. B716, 303 (2005), 
|hep-ph/0501042| . 

[27] M. Beneke, Y. Kiyo and K. SchuUer, Phys. Lett. B658, 222 r2008). CTw:0705.4518l 
[hep-ph] . 



31 



[28] M. Beneke, Y. Kiyo and A. A. Penin, Phys. Lett. B653, 53 (2007),[arXm07062733] 
[hep-ph] . 



[29] M. Beneke, Nucl. Phys. Proc. Suppl. 86 (2000) 547 |hep-ph/9910534] . 

[30] I. I. Y. Bigi, Y. L. Dokshitzer, V. A. Khoze, J. H. Kiihn and P. M. Zerwas, Phys. 
Lett. B 181, 157 (1986). 

[31] A. H. Hoang and C. J. Reisser, Phys. Rev. D 71, 074022 (2005), |hep-ph/0412258j. 

[32] M. Beneke, A. P. Chapovsky, A. Signer and G. Zanderighi, Phys. Rev. Lett. 93, 
011602 (2004), |hep-ph/0312331] . 

[33] M. Ben eke, A. P. C hapovsky, A. Signer and G. Zanderighi, Nud. Phys. B686, 205 
(2004), |hep-ph/0401002] . 

[34] M. Beneke, P. Falgari, C. Schwinn, A. Signer and G. Zanderighi, Nucl. Phys. B 792, 
89 (2008), irXiviOTOLOTTSj [hep-ph]. 

[35] S. Actis, M. Beneke, P. Falgari and C. Schwinn, Nucl. Phys. B 807, 1 (2009), 
larXiv:0807.0T02l [hep-ph]. 

[36] A. H. Hoang, C. J. Reisser and P. Ruiz-Femenia, Nucl. Phys. Proc. Suppl. 186, 403 
r2009). .arXiv:0810.2934l [hep-ph]. 

[37] A. Hoang, C. Reisser and P. Ruiz-Femenia, l arXiv: 1002,3223] [hep-ph]. 

[38] M. Beneke, in: Proceedings of the 8th International Symposium on Heavy Flavour 
Physics, Southampton, England, 25-29 July 1999, |h ep-ph/99 11490] . 

[39] M. Jezabek and J. H. Kiihn, Nucl. Phys. B 314, 1 (1989). 

[40] K. Melnikov and O. I. Yakovlev, Phys. Lett. B 324, 217 (1994), [hep-ph/ 93 02"31l] . 

[41] J. Alwall et ai, JHEP 0709, 028 f2007). larXiv:0706.2334l [hep-ph]. 

[42] M. Beneke and V. A. Smirnov, Nucl. Phys. B 522, 321 (1998) |hep-ph/971l'39ll . 

[43] J. A. M. Vermaseren, Comput. Phys. Commun. 83, 45 (1994). 

[44] D. Binosi, J. Collins, C. Kaufhold and L. Theussl, Comput. Phys. Commun. 180, 
1709 r2009). arXiv:0811.4113. [hep-ph]. 



[45] J. A. M. Vermaseren, |math-ph/0010025| . 



[46] R. Mertig, M. Bohm and A. Denner, Comput. Phys. Commun. 64, 345 (1991). 



32 



